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ABSTRACT 

Oscillations and propagating waves are commonly seen in high-resolution obser- 
vations of filament threads, i.e., the fine- structures of solar filaments/prominences. 
Since the temperature of prominences is typically of the order of 10"^ K, the promi- 
nence plasma is only partially ionized. In this paper, we study the effect of neutrals on 
the wave propagation in a filament thread modeled as a partially ionized homogeneous 
magnetic flux tube embedded in an homogeneous and fully ionized coronal plasma. 
Ohmic and ambipolar magnetic diffusion are considered in the basic resistive MHD 
equations. We numerically compute the eigenfrequencies of kink, slow, and Alfven 
linear MHD modes, and obtain analytical approximations in some cases. We find that 
the existence of propagating modes is constrained by the presence of critical values 
of the longitudinal wavenumber. In particular, the lower and upper frequency cut- 
offs of kink and Alfven waves owe their existence to magnetic diffusion parallel and 
perpendicular to magnetic field lines, respectively. The slow mode only has a lower 
frequency cut-off, which is caused by perpendicular magnetic diffusion and is signifi- 
cantly affected by the ionization degree. In addition, ion-neutral collisions is the most 
efficient damping mechanism for short wavelengths while ohmic diffusion dominates 
in the long- wavelength regime. 

Subject headings: Sun: oscillations — Sun: magnetic fields — Sun: corona — Sun: 
prominences 



1, INTRODUCTION 

Solar prominences and filaments are large-scale magnetic structures located in the solar corona. 
They are primarily composed of hydrogen (~90%) with a much smaller quantity of helium (~10%), 
their composition being similar to solar and cosmic abundances. The physical properties of their 
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plasma are very different from those typical of the coronal medium. The prominence temperature 
is of the order of lO'* K, i.e., a hundred times smaller than the coronal temperature, whereas the 
prominence density is around two orders of magnitude larger than that of the surrounding coronal 
plasma. For these reasons, the very existence of cold and dense filaments embedded in such an 
"extreme" coronal environment is not well-understood yet, but certainly the coronal magnetic field 
must play a crucial role in the formation, support against gravity, and thermal isolation of promi- 
nences. Due to their relatively cool temperature, the prominence plasma is only partially ionised. 
However, the exact ionization degree of prom inences is unknown and t he reported ratio of elec- 
tron density to neutral hydrogen density (e.g. JPatsourakos & Vialll2002h covers about two orders 
of magnitude (0.1 - 10). While the support of ions and electrons against gravity can be explained 
by the magnetic force, the support of neutral atoms seems more problematic because they are in- 
sensible to the magnetic field. Regarding th is problem, a plausible ex.planation for the support of 
neutrals is by means of collisions with ions (|Pecseli & Engvoldll2000l : [Gilbert et al.ll2002l . l2007|) . 
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clearly indicate that they are made by a myriad of thin (0".2 - 0".6) and long (5" - 20") thread- 
like structures piled up to form the filament body. These fine- structures, usually called fibrils or 
threads, seem to be orientated along magnetic field lines, and are observed in both active region 
and quiescent prominences. In observations from the Hinode spacec raft, threads usually appear as 



vertical structures in quiescent prominences (e.g. jBerger et al 



2008h. while ho rizontal threads are 



more commonly seen in active region prominences (e.g. JOkamoto et al.ll2007|) . However, this dis- 
tinction does n ot seem ve ry robust since the orientation of threads can sig nificantly vary within the 
same filament (jLinl 120041) . Theoretical models of filament threads (e.g.. iBallester & Priestill989l : 
Rempel et al.ll 19991) usually represent them as part of larger magnetic flux tubes whose footpoints 
are anchored in the solar photosphere. The large number of observations of oscillations and prop- 
agating waves in filament threads (s ee recent reviews by lOliver & Ballesteill2002l : lBallestei]|2006l : 
Banerjee et al.ll2007l : lEngvoldll2008b suggests that these phenomena are very frequent and ubiqui- 
tous in prominences. On the other hand, these oscillations are characterized by their quick damp- 
ing, the darn ping time being of the order of a few, typically four, periods (the reader is referred to 
Oliveiil2008L and references therein, for a review of this issue). 



Motivated by the observational evidence, a number of works studying filament thread oscilla- 
tions from a theoretical point of view have been performed within the last decade. By considering 
the magnetohydrodynamic (MHD) theory, early works studied the MHD eigenmodes supported by 
a filament thread modeled as a Cartesian slab, partially filled with prominence plasma, and embed- 
ded in the corona ( Joarder et al.lll997 : Dfaz et allboOl . 200 3 ). Later on, these works were extended 
by co nsidering a more representative cylindrical geometry (|Diaz et al. Ibooi iDvmova & RudermanI 
20051). It is worth mentioning that the collective oscillations of system s of threads hav e been also 



investigated in both Cartesian (IDfaz et al.ll2005l : iDiaz & Robertsll2006 i) and cylindrical dSoler et al 
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2009a() geometries. Focusing on t he damping of the oscillations, studies on t he effect of nonadia- 
batic mechanisms and mass flows (|Soler et alj|2008|) and resonant absorption (lArregui et al.ll2008 : 
Soler et alj|2009b|) have been performed. All these referred works neglected the presence of neu- 
trals and assumed a fully ionized hydrogen prominence plasma. Therefore, the effect of neutrals, 
and in particular that of ion-neutral collisions, on the propagation and damping of MHD waves in 
filament threads is not assessed yet and is the main motivation for the present investigation. 



There is an extensive literature regarding wave prop agation in a partia 



plasma in the context of laboratory plasrna physics (e.g.. 



Tanenbaum & Mintzerll 19621 : IWoods 



Watanabe 



196ilKulsrud & Pearcel[l969l : lwatts & Hannall2004l) . In as- 
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ally 
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trophysical plasmas, the typical frequency of MHD waves is much smaller than the coUisional 
frequencies between species. In such a case the one-fluid approach is usually adopted. One 
can find examples of works studying MH D waves in an unbounded, partially ioni zed, one-flui d 



plasma applied, e.g., to molecular clouds (lBalsaralll996l). to protoplanetarv disks (iDesch 



and to wave damp ing in the solar atmosphere (IDe Pontieu et al.l 12001 



Leake et al 



Khodachenko et al 



2005h. In t he con text of solar prominences, works by iForteza et al.l (|2007l |2008|) are 



2004), 
2004 



relevant. iForteza et al.l (l2007() derived the full set of MHD equations for a partially ionized, one- 
fluid plasma and applied them to stu dy the time darnping of linear, adiabatic waves in an unlimited 
prominence medium. Subsequently. IForteza et al.l (|2008|) extended their previous investigation to 
the nonadiabatic case by including thermal conduction by neutrals and electrons and radiative 
losses. Because of the effect of neutrals, in particular that of ion-neutral collisions, a generalized 
Ohm's law has to be considered, which causes some additional terms to appear in the resistive mag- 
netic induction equation in comparison to the fully ionized case. Among these additional terms, the 
dominant one in the linear regime is the so-called ambipolar magnetic diffusion, which enhances 
magnetic diff i ision across magnetic field lines. In this Paper, we apply the equations derived by 
Forteza et al.l (|2007|) to investigate the propagation of MHD waves in a partially ionized filament 
thread. Our filament thread model is composed of a homogeneous and infinite magnetic flux tube 
with prominence conditions surrounded by an unbounded and homogeneous coronal medium. The 
only prev ious investigation of wav es affected by ion-neutral collisions in a solar magnetic struc- 
ture is by Kumar & Roberts (|2003L but these authors considered the slab geometry and focused 
on surface waves in photospheric-like conditions. Hence, to our knowledge the present work is the 
first attempt to study MHD wave propagation in a partially ionized cylindrical flux tube. 

This paper is organized as follows: Section [2] contains a description of the model configura- 
tion and the basic equations. The results are presented and discussed in Section |3l Finally, our 
conclusions are given in Section |4l 
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2. MODEL AND METHOD 

2.1. Equilibrium Properties 

The model configuration considered in the present work is made of a homogeneous, isother- 
mal, and infinite plasma cylinder of radius a, representing a filament thread, embedded in a un- 
bounded coronal medium. We use cylindrical coordinates, namely r, ip, and z for the radial, az- 
imuthal, and longitudinal coordinates. In all the following expressions, a subscript indicates local 
values, while subscripts / and c explicitly denote filament (internal) and coronal (external) quan- 
tities, respectively. The magnetic field is taken homogeneous and orientated along the z-direction. 
Bo = BqZ, with 5() = 5 G everywhere. By adopting the on e-fluid approximation and considering 



only a hydrogen plasma (see details in iForteza et al.ll2007h . both the internal and external media 



are characterized by their plasma density, po, temperature, Tq, and number densities of neutrals, 
n„, ions, n,-, and electrons rig, with rig = rii. Thus, the gas pressure is = (Irii + n„)kBTo, where 
is Boltzmann's constant. The relative densities of neutrals, and ions, are given by 

= = (1) 

where we have neglected the contribution of electrons. We can now define a ionization fraction 
which gives us information about the plasma degree of ionization. 

This parameter h fx = 0.5 for a fully ionized plasma and jx = I for a neutral plasma. Any value 
of fx outside this range is physically meaningless. In the present work we take Tf = 8000 K, 
Pf = 5 X 10"'^ kg m~^, and Pc = p//200 = 2.5 x 10"'^ kg m~^. The coronal medium is assumed 
to be fully ionized, so jXc = 0.5. The filament ionization fraction, juf, is considered to be a free 
parameter, while the continuity of gas pressure selects the coronal temperature, 

T, = ^^Tf. (3) 

Pcfif 

Thus, the coronal temperature varies between Tc = 8 x 10^ K for /// = 0.5 and Tc = 1.6 x 10^ K 
for /ij = 1. MKS units are used throughout the paper. 



2.2. Basic Equations 



The detailed derivation of the basic MHD equations for a partially ionized hydrogen plasma 
in the one-fluid approach can be found in iForteza et al.l (|2007r) . After assuming small-amplitude. 
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adiabatic perturbations from the e quilibrium state, the basic equations are written in their linearized 
form (see Equations (18)-(21) of iForteza et alj|2007() . which in the notation adopted in the present 
work are 

dpi 



, Y7 



vi =0, 



>°o^ = -Vpi + -[(VxBi)xBo] 
dt fiQ 



dp 



dt ' dt ' 



~dr 



= V X (vi X Bo) + T/V^Bi + 7/aV X 



(V X Bi) X Bo 



XBn 



(4) 



(5) 



(6) 



(7) 



along with the condition V • Bi =0. In these equations Bi = [Br,B^,B^, \\ = (v^, v^^v^), pi, 
and pi are the magnetic field, velocity, density, and gas pressure perturbations, respectively, while 
cl = ypo/po and fiQ = 4n x 10"^ N A"^ are the sound speed squared and the vacuum magnetic 
permeability, with 7 = 5/3 the adiabatic index. Equations ©-([T]) are formally identical to their 
equivalent linear ideal MHD equations, with the exception of the induction equation (Equation (|7])) 
which contains additional terms due to the presence of neutrals and a nonzero resistivity. Parame- 
ters T] and ?7a are the coefficients of ohmic (also called Coulomb's) and ambipolar magnetic diffu- 
sion, respectively. In the present investigation, we neglect the Hall term in the induction equation 
because in a fully ionized plasma it may be important for frequencies larger than ~ lO'^ Hz, which 
are much larger than the observed frequencies of prominence oscillations. In a partially ionized 
plasma, the relative importance of the Hall effect gro ws with the density of neutrals, but in promi- 
nence conditions it can be still safely neglected (see iPandey & Wardlell2008l : iKrishan & Varghese 



20081). 



Equation ([7]) can be rewritten in a more convenient form to explain the role of the nonideal 
terms. For this purpose, we define the quantity S and express the ambipolar diffusivity as. 



(1 +^i)an 



IBn 



(8) 



where a,, is a friction coefficient and t/c is the so-called Cowling's coefficient of magnetic diffusion. 
Expressions for 77, 77c, and as functions of the equilibriu r n physical properties and the ionization 



fraction are found in the li terature (e.g., lBraginskiilll965l : iKhodachenko et al.l l2004l : iLeake et al 



2OO5I : iLeake & Arberll2006|) and are summarized in Appendix lAl Thus, Equation (|7]) becomes 
5Bi 



dt 



Vx(vi xBo)--Vxj||-— Vxj^-HVx(V;?i xBo). 

a (Tc 



(9) 



where cr = (/io?/) ' and (Tq = (MoVc) ' are the ohmic and Cowling conductivities, whereas jy and 
jx are the parallel and perpendicular components to the background magnetic field of the density 
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current perturbation, j, respectively, which are given by (|Arber et alj|2007h 

. a -Bo) Bo . BoXaxBo) 

J" = ^i^r ' = iBoP ' ^''^ 

with j = (V X Bi) /fio- Hence, we see that parameters rj and tjc correspond to the coefficients of 
magnetic diffusion parallel and perpendicular to magnetic field lines, while the term with the factor 
S is responsible for coupling the gas pressure gradient to the magnetic field evolution. For a fully 
ionized plasma, rjc = rj and S = 0, so that ambipolar diffusion is suppressed and the magnetic 
diffusion is isotropic. In such a case. Equation ^ becomes 

^ = Vx(vi xBo)-77VxVxBi, (11) 
ot 

which corresponds to the usual induction equation for a fully ionized, resistive plasma. However, 
in the presence of neutrals rjc ^ rj, meaning that the main effect of ion-neutral collisions is to cause 
magnetic diffusion to become a nonisotropic process. It is worth mentioning that rjc ^ r] even for 
a small relative density of neutrals. 

Since (p and z are ignorable coordinates, perturbed quantities are written proportional to 
exp (icot + irrnp - ik^z), where co is the frequency, and m and h are the azimuthal and longitudinal 
wavenumbers, respectively. After applying this Fourier-analysis to Equations (H])-©, one obtains 
the following system of eight equations which govern the propagation of linear MHD waves, 

/ Vr im \ 
lojpi = -po \v^ + — + — - ik.Vzj , (12) 



1 

iojVr = p\--^(iKB, + B'), (13) 

Po Bq ^ 

lojv^ = pi - [ik,B^ + —BA, (14) 

r Po Bq\ r I 

ik, 

iojv^ = —pu (15) 
Po 

ict>(j9i - c^pi) = 0, (16) 

(fTt^ lift lift \ / \ 

—Br + —B'^ + —BJ - 77c (k^Br - ik,B[) + ik,EBop[, (17) 

(„ 1 , 1 im , im \ / m o \ m 
K + -K - -^^v - —^r + TT^r -?7c(^z-5, + klB^yk,EBo-p,, (18) 

/ , 1 im \ /„„ 1 „, „ ., „, iky ^ , fn ^ \ 
iojB, = -Bo v' + -V,- + —v^ + 77c 5 ' + -Bi - —B, + ik,B' + - k,-B^ 
\ r r I y ^ f ^ f2 ' f. r I 

+S5o(k + Vi-^Pi), (19) 
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where the prime denotes derivative with respect to r and v\ = 5q///oPo is the Alfven speed squared. 
Equations (fT2l)-(fT9l) form an eigenvalue problem which we numeric ally solve by means of the 
PDE2D code (ISewellll2005|) based on finite elements (see, for example. iTerradas et alJIlOOSL for an 
explanation of the method). The numerical integration of Equations (fT^ -(fT91) is performed from 
the cylinder axis, r = 0, to the finite edge of the numerical domain, r = r^ax- The evanescent 
condition is imposed in the coronal medium, so all perturbations vanish at r = r^ax- We impose 
the evanescent condition because we focus our study on trapped, nonleaky modes supported by 
the filament thread body. It is unlikely that the observed thread oscillations correspond to leaky 
mode s, which are high-frequency propagating waves in the coronal medium (see for example lCally 



1986|) . while an interpretation in terms of trapped solutions is more realistic and consistent with 



the observed periods. Therefore, the edge of the numerical domain has been located far enough 
from the filament thread to obtain a good convergence of the solution and to avoid numerical 
errors (typically, we consider Vax > 100a). The boundary conditions at r = are imposed by 
symmetry arguments. We assume a fixed, real, and positive k^, so the numerical solution provides 
with a complex oscillatory frequency, co, as well as the eigenfunctions of perturbations. Wave 
solutions appear in pairs, a)+ = +ojr + ico\. The solution a>+ corresponds to a wave propagating 
towards the positive z-direction, whereas the contrary stands for Since there are no flows in the 
equilibrium configuration, both solutions have the same properties, and so we restrict ourselves to 
solutions with % (a») > 0. The oscillatory period, P, the damping time, td, and the damping ratio, 
td//', are related to the frequency as follows 

n 2;r 1 Td 1 a;R 

P= — , Td = — , — = . (20) 

a»R a»i F in (x)\ 



2.3. Dimensional Analysis of the Induction Equation 



Forteza et al.l (120071 see their Figure 6) evaluated the relative importance of the nonideal 
terms of the induction equation for typical prominence conditions. For the range of parameters 
considered in their paper, these authors concluded that the effect of the pressure-related term, i.e., 
that with the parameter S, is negligible, while the term with Cowling's diffusivity is the dominant 
one. Here we quantify the importance of ohmic and Cowling's diffusion by means of a dimensional 
analysis. Equation ^ indicates that magnetic diffusion in the parallel and perpendicular directions 
is governed by the terms with ohmic diffusivity and Cowling's diffusivity, respectively. We thus 
define the parallel and perpendicular magnetic Reynolds numbers as follows, 

r]L T]cL 

where L^j and L^,^ are typical length-scales related to the terms of the induction equation with rj and 
rjc, respectively, whereas L = min (^Lrj,L^^y and C/q is a typical plasma velocity. In the absence of 



Rm\\ - Rm± - (21) 
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neutrals, L^, = L^^^ = L, so Rm\\ = R,n± = because magnetic diffusion is isotropic. However, 
since this mechanism is nonisotropic in a partially ionized plasma, both the longitudinal and per- 
pendicular directions have different typical length-scales. An examination of Equations (fTTI) . (fTSi) . 
and (fT9l) reveals that longitudinal derivatives only appear in the terms with rjc, whereas those with 
i] only contain radial and azimuthal derivatives. Hence, Ljj is a typical length-scale perpendicular 
to magnetic field lines and can be related to the filament thread radius, ~ a. On the other hand, 
Lrj^ corresponds to a typical length-scale along the magnetic field direction and can be expressed 
in terms of the longitudinal wavelength, Lrj^ ~ or, equivalently, in terms of the longitudinal 
wavenumber, Ljj^ ~ 2nk~^. In addition, we relate the velocity- scale to the filament thread sound 
speed, Uq ~ Cs/. Therefore, the Reynolds numbers (Equation (|2TI) ') are rewritten as follows, 

Rm\\ = J Rml. = -p^- (22) 

We see thati?„,|| is independent of the longitudinal wavenumber, while 7?m± is inversely proportional 
to kl. This suggests that the relative importance of Cowling's diffusion increases with h. To 
perform a simple calculation we take k^a ~ 1, consider a typical value a = 100 km, and assume 
fif = 0.8. We obtain R^w « 7 x 10^ and Rmi. ~ 4x 10^. For the considered parameters, the dominant 
nonideal mechanism is Cowling's diffusion, while ohmic diffusion has a minor role. Furthermore, 
it is possible to obtain an estimation of the wavenumber for which both ohmic and Cowling's 
diffusion have the same importance by equaling the parallel and perpendicular magnetic Reynolds 
numbers, i.e., R^w = R,n±. By this procedure we obtain, 

k^a ~ 2n ^ — , (23) 

which gives k^ a ~ 2.8 x 10 " ^ for the same parameters assumed before. These estimations are 



consistent with Forteza et al.l (|2007|) . since they considered quite a large value for the longitudinal 



wavenumber and therefore ohmic diffusion is not relevant in their case. We have to bear in mind 



that the observed width of filament threads is in the range 0".2 - 0" .6 (|Linll2004|) . and therefore 
a ranges from 75 km to 375 km, approximately. On the other hand, the detected wave lengths of 
prominence oscillations are between 5 x lO"* km and 10^ km ([Oliver & Ballesteiil2002|) . One can 



combine both quantities (the wavelength and the thread width) into the dimensionless quantity k^a 
and compute its upper and lower limits. So, for a = 75 km, 5 x 10"^ < k^a < 9 x 10"^, while for 
a = 375 km, 2x 10"^ < ha < 4x 10 ^ Thus, taking into account both intervals and considering that 
thinner threads than those resolved by present-day telescopes might exist, the relevant range of ha 
of prominence oscillations covers two orders of magnitude and corresponds to lO""* < k^a < 10"^ 
This range of k^a contains all realistic values of the wavelength and the thread width. According 
to Equation (|23]) . both ohmic and Cowling's diffusion could be important in such a range of k^a. 
This will be verified through the numerical computations in the next Sections. 
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3. RESULTS 



3.1. Alfven Waves 



We start our investigation by studying Alfve n waves. Some works have investigated Alfven 



wave propagation in partially ionized plasmas (e.g.JWatanabe 



Watts & Hannall2004 IPandev & Wardld l2008l : Forteza et al 



1961bl : lTanenbaum & Mintzej 19621 : 



20081) . but to our knowledge a de- 



tailed investigation in cylindrical geometry applied to solar plasmas remains to be done. In general 
Alfven modes are coupled to magnetoacoustic modes except for m = 0. Thus, here we assume no 
azimuthal dependence in order to study Alfven waves separately. We see that by setting m = 
Equations (fT4l) and (fTSl) are decoupled from the rest, 



iiov^ = -ih—B^, 
Bq 



iojB,, 



-ik,BoV^ + r]\B'' + -B' -\bJ- r]cKB^. 



(24) 



(25) 



r ^ 

These equations only involve the perturbations and B^, so Alfven modes are purely torsional, 
incompressible waves for m = 0. Now, we use Equation (l24l) to express as a function of B^, and 
after substituting it in Equation (|25l) a single equation for B^ is obtained. 



tjB; + Tj-^B'^ + 



B^ = 0. 



Equation (|26]) can be rewritten as follows, 

r^B; + rB'^ + (my-l)B^ = 0, 



with 



m. 



(26) 



(27) 



(28) 



The same equation stands for v^. So, the Alfven wave is governed by a Bessel equation of order 
1, where plays the role of the r adial wavenumber an d Fa is equivalent to the modified Alfven 
speed defined in Equation (29) of iForteza et al.l (|2008[ ). We see then that in the resistive case, 
Alfven modes ar e not strictly confined to rnagnetic surfaces and becorn e "global" eigenmodes of 
the flux tube (see lFerraro & Plumptonlll96ll : ISylll984l : ICopil et al.ll2008|) . The general solution of 
Equation (|T7]) for regular perturbations at r = and vanishing at infinity is 



B^(r) 



Ai7i(mA/r) 



if 
if 



r < a, 
r > a. 



(29) 



Al and A2 being complex constants. 7i and //[^^ are the usual Bessel func tion and Hankel function 
of the first kind, respectively, of order 1 (|Abramowitz & Stegunlll972h . In order to obtain the 
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dispersion re lation, we need to impose boundary conditions at r = a. iTomimural (|1990|) (see also 
Woodslll962|) provide appropriate boundary conditions for our case, 



(30) 



where [[X]] = Xc - Xf stands for the jump of the quantity X at r = a. is set continuous across 
the boundary because we assume no surface currents on the interface between the thread and the 
corona. On the contrary, the jump of the derivative of B^ is determined by the difference between 
the internal and external ohmic diffusivities. Applying these boundary conditions and after some 
algebra, the following dispersion relation is obtained. 



(niAca) J[(mAfa) 
rifMAf- 



nf - Tic 



Hf\mAca) 



(31) 



JiiiTiAfa) a 

Next, we solve Equation (|3T1) and plot in Figure [IJa) the Alfven wave phase speed as a function 
of k^a for several ionization degrees. As one can see, the Alfven mode only propagates between 
two critical wavenumbers, the first of them being independent from the ionization degree. At these 
critical wavenumbers, the real part of the Alfven frequency vanishes. Note that these critical values 
are far from the range of k^a relevant for prominence oscillations. On the other hand. Figure (HZ?) 
displays the ratio td//' as a function of k^a. tu/P is independent from /ij for small k^a, while it is 
significantly afl^ected by the ionization degree for large ha. The maximum of Td/P occurs within 
or close to the relevant range of k^a. Figure[2lallows a better understanding of the behavior of td /P- 
Here, we display a comparison of the actual value of td/P with that obtained by neglecting one 
of the two possible damping mechanisms, i.e., ohmic diffusion (by setting 77 = 0) or ion-neutral 
collisions (by setting rjc = rf). As expected, the presence of neutrals has an important effect on 
the damping time for large k^a, whereas ohmic diffusion dominates for small ha. As estimated 
by Equation (|23]) . the transition between both regimes takes place within the relevant range of k^a, 
and the estimated transitional wavenumber is close to the actual value. 

To perform a more in-depth study of Alfven wave propagation, let us consider now the case 
in which ohmic diffusion can be neglected, i.e., k^a much larger than the transitional one (Equa- 
tion (I23])). So, by taking 77 = in Equation (|26l) one obtains that Alfven waves verify. 



y [klv\ - 0?) - lick, 

CO ^ ' 



B^ = 0, 



(32) 



which for an arbitrary B^ gives us the dispersion relation for Alfven modes in the regime dominated 
by Cowling's diffusion. 



o) - irjck CO 



k'A = 0. 



(33) 



Equation (1331) is formally identical to Equation (49) of IPandey & Wardld (|2008|) . Its exact solution 
gives us a complex frequency, 

K 
2 



CO 



^44 



(34) 
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Fig. 1. — Results corresponding to the Alfven mode, (a) Phase speed as a function of ha in units 
of the internal Alfven speed, (b) Ratio of the damping time to the period as a function of ha. 
In both panels, different linestyles represent different ionization degrees: jlf = 0.5 (dotted line), 
/// = 0.6 (dashed line), /// = 0.8 (solid line), and fif = 0.95 (dash-dotted line). Symbols are the 
approximate solution given by Equation (l36l) for /i/ = 0.8. The shaded zone corresponds to the 
range of typically observed wavelengths of prominence oscillations. 
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Fig. 2. — Ratio of the damping time to the period as a function of ha corresponding to the Alfven 
mode with \Xf = 0.8. The solid line is the complete solution considering all terms in the induction 
equation, i.e., the solid line in Figure[IlZ?). The dotted and dashed lines are the results obtained by 
neglecting ohmic diffusion (77 = 0) or ion-neutral collisions (rjc = rf), respectively. The vertical 
dot-dashed line is the analytically estimated transitional wavenumber between both regimes of 
dominance (Equation (l23l)). 
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suit that the Alfven wave is damped in a 


partially ionized plasma (see, 


e.g.. 


Haerendel 


1992; 


De Pontieu & Haerendel 


1998 




De Pontieu 


1999 




De Pontieu et al. 


2001 




Leake et al. 


200 


5, for 



studies on the damping of Alfven waves in the solar chromosphere). Obviously, if t/c = we 
recover the ideal, undamped Alfven mode, oj = +k^VA- From Equation (|34l) we can also see that 
the real part of the Alfven frequency vanishes for a critical value of h. 



K = 



2va 



(35) 



This critical wavenumber is equivalent to that given by Equation (38) of iForteza et al.l (|2008r) by 
considering parallel propagation. For > the Alfven wave becomes a nonpropagating, purely 
damped disturbance. This critical wavenumber corresponds to the largest critical wavenumber 
obtained in Figure [T] 

To obtain an expression for the smaller critical wavenumber, we have to consider the case 
T] 0. To make some analytical progress we start by plotting the perturbation corresponding 
to the radially fundamental Alfven mode (Figure [3]). Although the eigenfunction is not strictly 
confined within the cylinder, we note that the amplitude of for r > a is much smaller than 
that within the thread. This is caused by the large contrasts of the internal and external densities, 
PflPc = 200, and ohmic diffusivities, rff/ric ~ 3 x 10^. On the basis of this evidence, it is possible 
to give an analytical approximation of the Alfven frequency in the general 77 ^ case by neglecting 
the influence of the corona and setting 5^ « at r = a. This implies that mp^fU ~ ji, with 71 the 
first zero of the Bessel function 7i. In this approximation, the Alfven mode frequency is 



OJ 



4 a/4^a^' 



T]ck^ + 7]^ 



a" 



(36) 



This expression is consistent with previous results since for 77 = it reduces to Equation (|34l) . As 
before, we can obtain an expression for the critical k^ at which the real part of the Alfven frequency 
vanishes. 



(37) 



We see that two fc^ are now possible. That given by the + sign, namely k^'^, is a correction of 



the previously described critical wavenumber for 77 = (Equation (pSl) ). On the other hand, that 
given by the - sign, namely k^f, is a new critical value which arises because of the combination 
of two effects, namely the nonzero value of j] and the presence of the factor given by the 



geometry. Since iForteza et al.l (|2008[) considered an infinite medium, this geometry -related critical 



wavenumber is absent in their investigation. Considering that rjcrifja^vl 1, a first-order Taylor 
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expansion of Equation (|37]) gives 



T/c \ 2v\ I rjc 2va 

One can see that k';~ does not depend on t/c and so it is not affected by the plasma ionization 
degree. Furthermore, since k';' < k^^*, the Alfven wave only exists as a propagating mode for 
k'- <k,< k'^. 

Z ^4 



1 .2 




Fig. 3. — Modulus in arbitrary units of the eigenfunction as a function of r/a corresponding to 
the radially fundamental Alfven mode with m = 0, jxf = 0.8, and k^a = 10"^. The vertical dotted 
line corresponds to the filament thread edge. The small panel shows a close-up of the eigenfunction 
close to the boundary between the thread and the corona. 

For fif = 0.8 the approximation given by Equation (|36l) is compared with the complete solu- 
tion by means of symbols in FigurefU We see that there is a good agreement between both solutions 
meaning that both the approximation to the frequency (Equation (|36l) ) and the expressions for the 
critical wavenumbers (Equation (|38l) ) are consistent with the actual result. The Alfven wave phase 
speed is identical to va in a fully ionized, ideal plasma. Therefore, we can perform an analogy and 
define an effective internal Alfven speed from the real part of Equation (|36l) . 

We see that except when k^ is close to one of the critical wavenumbers, the phase speed is almost 
equal to va/. 
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3.2. Kink Mode 

Here, we turn our attention to the kink mode. We fix the azimuthal wavenumber to m = 1 
since we are interested in modes that are able to perturb the cylinder axis and are reasonable can- 



didate s to produce the observed transverse oscillations of prominence threads (e.g. JOkamoto et al. 



20071) . We focus on the radially fundamental kink mode. It is not possible to give an analytical 
dispersion relation for this wave when both ohmic and ambipolar magnetic diffusion are consid- 
ered in the induction equation. For this reason, we numerically solve the full eigenvalue problem 
(Equations (fT2l) -(fT9l)) and obtain the frequency with the PDE2D code. 

Figure|3a) displays the kink mode phase speed as a function of k^a. We see that the propaga- 
tion of the kink wave is also constrained by the existence of two critical wavenumbers, that turn out 
to be the same found in the case of the Alfven wave, see Section [3TT1 This behavior can be easily 
understood by c onsidering the first-order asymptotic approximations of the kink mode frequency 



in the ideal case (|Edwin & Robertslll983h . 



^ _ I Ckk, for k,a<^l, ^^^^ 
\ v^fk^ for fc^a » 1, 

where Ck = yJrT^J^^^f '•^^ so-called kink speed. Note that for pf » pc, q ~ V2va/. So, 
both in the short- and long-wavelength regimes the kink mode frequency directly depends on the 
Alfven speed. Replacing the ideal Alfven speed by its effective value, va/, from Equation (|39l ), one 
sees that the kink mode frequency vanishes at the same critical wavenumbers as the Alfven wave. 
Drawing our attention within the relevant range of k^a, we obtain co/k^ « Ck, so nonideal effects do 
not seem relevant to wave propagation for the observed wavelengths. Next, the dependence of the 
ratio of the damping time to the period with k^a is shown in Figure SfZ?). This result is very similar 
to that obtained for the Alfven mode, hence no additional comments are needed in this case. 

Some analytical progress can be performed by adopting the 13 = case and neglecting ohmic 
diffusion. In such a situation, from Equation ^ one can write the magnetic force term as, 

(VxB,)xBo = ^Vi, (41) 

which can be used to rewrite Equation ([7]) in a compact form as follows. 



^Vx(vjXBo), (42) 



dt v\ 

with defined in Equation (|28l) . Equation ((42)) is formally identical to the ideal induction equation 
with the extra factor T\lv\ in the convective term. Note that ^\lv\ = 1 for rjc = 0. Thus, it is now 
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Fig. 4. — Results corresponding to the kink mode, (a) Phase speed as a function of k,a in units 
of the kink speed, {b) Ratio of the damping time to the period as a function of k^a. The different 
linestyles and the shaded zone have the same meaning as in Figure [U Symbols are the approximate 
solution given by solving Equation (1431) for jlf = 0.8. 

straight-forward to obtain a dispersion rel ation by following the same procedure as in the ideal 
case (see details in lEdwin&Robertsll 19831) . 



Pctnf [oj - kj^^) — = pfM, [co - kj:^^) (43) 

where K„, is the modified Bessel function and the quantities rrif and m,. are given by 

m) = '—-^ 1^, ml = '—-^ (44) 

^ A/ A Ac 

Note that Equation (l43l) applies to any value of m. Equation (|43l) has been solved for m = 1 and 
jlf = 0.8, the corresponding results being the symbols in Figure IH This approximate analytical 
solution agrees well with the numerical solution in the range of La dominated by Cowling's diffu- 
sion. For small kM, the numerical and analytical solutions do not agree because the effect of ohmic 
diffusion is missed by the analytical approximation. 

Next, we check here the efficiency of the nonideal terms of the induction equation on the kink 
mode damping. Figure [5] shows a comparison of the value of td/P of the kink wave with that 
obtained by neglecting ohmic diffusion (77 = 0) or ion-neutral collisions (77c = tj and S = 0). Only 
the solutions corresponding to /i/ = 0.8 are displayed. The result for the kink mode is very similar 
to that obtained for the Alfven mode (see Figure O: ohmic diffusion dominates for small k^a and 
ion-neutral collisions are dominant for large k^a. The transition between both behaviors takes place 
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within the relevant range of k^a, where ty)IP ^ 4. Therefore, neither ohmic diffusion nor ion- 
neutral collisions can provide damping times compatible with those observed. Only for an almost 
neutral plasma (fif > 0.95) and large L, one can obtain a Tu/P consistent with the observations. 
In comparison with the damping mechanisms studied in previous works, the efficiency of ohmic 
diffusion and io n-neutral collisions is much snialler th an that of resonant absorption for the kink 
mode damping dArregui et al.ll2008l: Isoler et al.ll2009bl) . 




Fig. 5. — Same as Figure[2]but for the kink mode. 



3.3. Slow Mode 

The last wave mode that we study in this investigation is the slow mode. Again we numerically 
obtain the frequency of the radially fundamental mode. We consider m = 1 but we note that the 
slow mode behavior is weakly affected by the value of the azimuthal wavenumber. Figure [6ta) 
shows the slow mode phase speed as a function of k^a. The slow mode behavior is also affected 
by the presence of a critical wavenumber which is highly dependent on the ionization degree. The 
slow mode is totally damped for k^ smaller than the critical value. We see that for large enough fif, 
the critical h falls inside or is larger than the observed values. This result has relevant implications 
from the observational point of view, since it suggests that the slow wave might not propagate 
in realistic, thin filament threads. For ha larger than the critical value, the slow mode phase 
speed is close to the internal tube (or cusp) speed, ct/ = Cs/Va/^/ ■yj'^lf + ^a/- other hand. 

Figure [6tZ?) displays Td/P again as a function of ha, where one can see that the damping time 
decreases dramatically as k^a approaches the critical value. Also, we see that the damping time 
grows rapidly when /// ^ 0.5, i.e., the dependence on the ionization degree is more sensitive to 
the ionization degree for an almost fully ionized plasma than for a weakly ionized plasma. 
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10'° 10'* 10'^ 10° 10^ 10'° 10'* 10'^ 10° 10^ 

k^a k^a 



Fig. 6. — Results corresponding to the slow mode with m = I. (a) Phase speed as a function of k^a 
in units of the internal cusp speed speed, (b) Ratio of the damping time to the period as a function 
of La. The different linestyles and the shaded zone have the same meaning as in Figure[T] Symbols 
correspond to the analytical approximate solution for jlf = 0.8 (Equation (|46|)'). 



Such as for the Alfven and kink waves, it is possible to perform some simple analytical calcu- 
lations in order to get a better understanding of the slow mode behavior. Here, it is crucial to take 
into account the common result regarding the i deal slow mode pr opagation in both the Ca rtesian 
slab JEdwin & Robertslll982l : Isoler et all r2007h and the cylinder fedwin & Robertslll983l) . It is 
well-known that the slow mode in a jS < 1, homogeneous medium is mainly polarized along the 
magnetic field direction. So, if one considers a magnetic structure, say, a slab or a cylinder, the 
slow mode is almost insensible to the perpendicular geometry to magnetic field lines and is mostly 
governed by the physical conditions internal to the structure. In such a case, the slow mode fre- 
quency is well approximated by solving the dispersion relation corresponding to a homogeneous 
medium with the physical conditions of the internal structure, and selecting an appropriate value 
for the perpendicular wavenumber to the magnetic field, namely kj^. This kj^ contains the effect of 
the geometry of the magnetic structure. Thus, to apply this technique we co nsider the dispersion 
relation for a partially ionized, homogeneous, and infinite plasma derived by iForteza et al.l (|2007l 
Equation (24)), which in our notation is 

co'^ - i [k] + k]) rico? - [c] + vy [k] + k]) a? 

+ic] (kl + kl) [[k] + kl) rjc - klEpovl] oj + (kl+ kl) klc]v\ = 0. (45) 

Note that in Forteza et al.l J2007|) the magnetic field is orientated along the x-direction, so their L 
and k^ correspond to our k^ and k^, respectively. We are interested in approximating the slow mode 
frequency near the critical k^, where oj is small and tends to zero. In such a situation, it seems 
reasonable to neglect the terms with or' and oJ^ in Equation (l45l) . Also, we neglect the term with 
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S. So, the dispersion relation becomes a second order polynomial for the frequency that can be 
solved exactly. 

We check the validity of this approximation by taking t/c = 0, so a» « icjk^, which is consistent 
with the slow mode frequency in the ideal case. The critical wavenumber, k^^'\ can be derived by 
setting the real part of Equation (|46|) equal to zero, 

K' = — ^\h^-kl. (47) 

In principle, we obtain two critical values from Equation (l47l) . However, that given by the + sign 
is not consistent neither with the small-frequency approximation nor with our numerical results 
because it corresponds to a very large value of k^. On the contrary, that given by the - sign is the 
critical value that we are looking for, and for kj_Cjri^/v^^ 1 it can be approximated by 

« (48) 

Hence, the critical wavenumber is both determined by the geometry, through k^, and by the ion- 
ization degree, through t/c- As a consequence, the slow mode has no critical wavenumber in the 



unbounded medium, studied by lForteza et alj (|2008|) 



To perform a comparison between the analytical approximation (Equation (|46l )) and the nu- 
merical results, we have to provide a value for k^^. The dominant velocity perturbation for the 
slow mode is v^. Again, we consider the well-known behavior of perturbations in the ideal case 
and assume that, for m = 1, Vj. is approximately proportional to the Bessel function 7i. Let us 
also assume that the corona has a negligible effect and the slow mode is fully confined within the 
filament thread. In such a situation, we can approximate k\ ^ j\la^ + ni^/a^, where the term j\/a^ 
is the radial contribution to the perpendicular wavenumber, while in^/a^ represents the azimuthal 
contribution. We take fif = 0.8 and compute the approximate slow mode frequency from Equa- 
tion (|46l ). This result is plotted by means of symbols in Figure [6l We see that the approximation 
is reasonably good for k^a < 1 and near the critical L, so the approximation describes well the 
presence of a critical wavenumber for the slow mode. However, it diverges from the numerical 
result for ha > 1 because the assumptions considered for deriving Equation (|46l) are not satisfied 
for large k^a. 

It is possible to know the relative position of the slow mode critical wavenumber (fc^* from 
Equation (|48l) ') with respect to the Alfven mode lower critical wavenumber (/r^" from Equation (IBS])). 

k" « ^^V-. (49) 

T] Va 
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For a partially ionized plasma rjc » rj, so we expect k^J » k^~. However, in the fully ionized case, 
rjc = J] but cj < Va, therefore kf/ < fc^". This is in agreement with our results, e.g., compare the 
lower critical wavenumbers in Figures \lla) and Sfa), on one hand, and Figure [6ta), on the other 
hand. 

Finally, the particular effect of ohmic diffusion and ion-neutral collisions on the the slow 
mode damping is assessed in Figure U\ We obtain that slow mode damping is entirely governed 
by ion-neutral collisions. The ratio Tu/P of the slow mode achieves very small values close to 
the position of the critical wavenumber. This is because the real part of the frequency tends to 
zero and so the period tends to infinity. It is worth noting that ohmic diffusion never becomes 
important for the slow mode damping because the presence of the critical wavenumber does not 
allow the slow mode to propagate for small values of k^a. In addition. Figure [6tZ?) indicates that in 
the case of strong ionization, the ion-neutral collision damping is less effi cient than the dam ping by 
nonadiabatic effects (thermal conduction and radiative losses) studied in lSoler et al.l (|2008|) . Only 
when juf is large, ion-neutral collisions become more efficient than nonadiabatic mechanisms for 
the slow mode damping. 




k,a 



Fig. 7. — Same as Figure[2]but for the slow mode. The dotted and solid lines are superimposed. 



4. CONCLUSION 



In this paper, we have studied the propagation of MHD waves in a partially ionized filament 
thread. Contrary to the fully ionized, ideal case, we have found that wave propagation is con- 
strained by the presence of critical values of the longitudinal wavenumber. This result is of special 
relevance in the case of the slow magnetoacoustic mode, since the critical wavenumber lies within 
the observed range of wavelengths of prominence oscillations. This might cause the slow wave 
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not to propagate in thin filament threads and, therefore, it might not be observationally detectable. 
Turning to the wave damping, we have obtained that both Alfven waves and magnetoacoustic 
waves are damped by ion-neutral collisions and ohmic diffusion. In the case of Alfven and kink 
waves, ion-neutral collisions (by means of Cowling's diffusion) dominate for short wavelengths, 
while ohmic diffusion is more important for large wavelengths. On the contrary, the slow mode is 
totally governed by Cowling's diffusion. With the exception of the slow mode, theses mechanisms 
cannot provide damping times compatible with those observed when typical values of the wave- 
length are considered. In the case of both Alfven and kink waves, one has to consider very short 
wavelengths and almost neutral plasmas to obtain realistic values of the damping time. 

One must be aware that in this work we have considered a plasma composed only by hydro- 
gen. So, it is worth to estimate the eff'ect of the presence of helium on the results. Helium remains 
mostly neutral for typical prominence temperatures, and therefore its presence would increase the 
value of p.f and enhance the efi'ect of ion-neutral collisions. As a consequence, the transitional k^a 
at which ion-neutral collisions become important might be shifted toward smaller values, mean- 
ing that Cowling's diffusion may dominate for longer wavelengths and/or thinner threads than in 
a pure hydrogen plasma. Also, the damping time of the three wave modes could be reduced by 
the presence of helium. On the other hand, since helium only corresponds to 10% of the promi- 
nence material approximately, it could have a minor influence, although its precise effect should 
be assessed in detail by future investigations. 

This paper is a first step in the study of MHD waves in partially ionized filament threads. 
There are several issues of special interest that may be broached in future works. Some of them 
are, for example, to assess the effect of neutrals on the resonant absorption of kink waves in the 
Alfven continuum, to study the spatial damping of waves, and to investigate the time-dependent, 
initial value problem. 
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funds, and the Conselleria d'Economia, Hisenda i Innovacio of the CAIB under Grants No. AYA2006- 
07637 and PCTIB-2005GC3-03. RS thanks the CAIB for a fellowship. 



A. EXPRESSIONS FOR THE PARAMETERS USED IN THIS PAPER 

We give here expressions for the quantities rj, tjq, and Qr„ that appear in the text. All quantities 
are expressed in MKS units. Ohmic and Cowling's difFusivities are given by Coulomb's, cr, and 
Cowling's, cTc conductivities, respectively, 

n = — > nc = , (Al) 
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with 

where and e are the electron mass and charge, and is the friction coefficient, given by 

Qr„ = meHeV^^ + m^n/y^^. (A3) 

There V-, v^„, and are the effective electron-ion, electron-neutral, and ion-neutral coUisional 
frequencies, respectively, defined as 

Vei = ; ^eh Ven = \ "^en, ^in = \ Vin, (A4) 

m,- + me rnn + me mn + m, 



with 



_f,nihZ?- _ ShTo _ UbTo 

Vei = 3.7 X 10 ^en — \ ^enj ^/n = ^In -\ ^/n; (A5) 

where 2,„ = 5x10"!^ and2,„ = 10" are the ion-neutral and electron-neutral cross-sections, 
respectively, A is de Coulomb logarithm, and 

m,„ = , men = • (A6) 

mi + m„ me + mn 

For a hydrogen plasma, Z = 1 and m, ~ m„, so a simplified expression for the friction coefficient 
an can be provided by neglecting the contribution of electron-neutral collisions. 



an = -^n(i- ^n) A ^in- (A7) 

2 m„ ^ nmi 

In a fully ionized plasma, ^„ = 0, hence crc = <t and Qr„ = 0. Also note that in a fully neutral 
plasma, ^„ = 1 and n, = = 0, so cr = crc = Qr„ = 0. 
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